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Abstract 

We present a QCD sum rule calculation of the Isgur-Wise form-factor t\(v ■ 
v') and ti(v ■ v') for the semileptonic decays B — > Di(2420)£^ and B — > 
Z?2(2460)£p in the framework of heavy quark effective theory. These two 
universal functions, associated with the matching of the weak currents in 
QCD onto those in the effective theory, appear at the order 1/mg in the 
heavy quark expansion of meson weak decay form factors. 
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I. INTRODUCTION 



The heavy quark effective theory (HQET) has important consequences on the 
spectroscopy and the decay matrix elements of the hadrons containing a single heavy 
quark Q. In the infinite mass limit, the spin and parity of the heavy quark and that of 
the light degrees of freedom are separately conserved. This allows that the hadronic states 
can be classified in degenerate doublets by the total angular momentum j and the angular 
momentum of the light degrees of freedom ji- in the case of qQ mesons, coupling ji with 
the spin of heavy quark sq = 1/2 yields a doublet with total spin j = ±1/2. The ground 
state mesons with jf = | are the doublet (D,D*) for Q = c and (B,B*) for Q = b. The 
excited heavy mesons with jf = 1+/2 and 3 + /2 can be classified in two doublets of spin 
symmetry (0 + ,l + ) and (l + ,2 + ), respectively. The charmed mesons with jf = 3 + /2 have 
been experimentally observed. For q = u,d, these mesons are denoted as the -Di(2420) and 
the 1^2(2460) mesons. While the jf = l + /2 doublet (D' Q , D[) has not been observed yet. 

One of the most important applications of the heavy quark symmetries has been the 
study of semileptonic transitions between two heavy hadrons. The hadronic matrix ele- 
ments of weak currents between members of the doublets identified by je and jgi can be 
expressed in terms of universal form factors which are funtions of the dot-product, y = v-v', 
of the initial and final hadron four-velocities. A well-known result is that the semileptonic 
decays B — ► Dlv and B — > D*W, in the uiq — > 00 limit, can be described in terms of a 
single universal function, the Isgur-Wise function £(y). In the case of transitions between 
members belonging to different heavy quark multiplets, additional form factors needs to be 
introduced. For example, for the B semileptonic decay into excited charmed mesons (D' Q , 
D'i) and (Di, D2) the weak matrix elements can be expressed in terms of two independent 
functions, r(y) and £($/) [§|,|J, respectively, in the limit uiq — > 00. 

There are Aqcd/^q corrections to the weak matrix elements parametrized by form 
factors at the niq — > 00 limit. The Aqcd/^q corrections to the leading term can be 
analyzed in a systematical way in HQET. The matrix elements receive contributions from 
higher-dimensional operators in the effective currents and in the effective Lagrangian. They 
result in a set of new Isgur-Wise functions. 

The universal functions must be estimated in some nonperturbative approaches. A vi- 
able approache is the QCD sum rules formulated in the framework of HQET. This method 
allows to relate hadronic observables to QCD parameters via the operator product ex- 
pansion (OPE) of the correlator. A fruitful application of QCD sum rules has been the 
determination of the Isgur-Wise functions parameterizing the B — ► semileptonic tran- 
sitions up to the AQ CD /mQ corrections P-|7|]. 

The semileptonic B decays into excited charmed mesons have attracted attention in 
recent years. From the phenomenological point of view, the B — > D** transitions are inter- 
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esting (here D** denotes the generic L = 1 charmed state), since in principle these decay 
modes may account for a sizeable fraction of the inclusive semileptonic -B-decay rate. At 
leading order in the 1/tuq expansion, the two independent universal form factors, r(y) and 
C(y), that parametrize the transitions B — > (D\,D?£)lT> and B — > (D' Q , D[)£i>, respectively, 
have been calculated with QCD sum rules [|]|§. Moreover, perturbative corrections to 
0(a s ) have been included in the QCD sum rule for ((y) in fTOfl . The other approaches 



include different quark models Jl I {I7J and relativistic Bethe-Salpeter equations [IS|. A 



model independent analysis has been carried out in [|J for the Aqcd/^q corrections. At 
the order 1/rriQ, the corrections for matrix elements of B — > D** include contributions 
from higher-dimensional operators in the effective currents and in the effective Lagrangian. 
For the semileptonic transitions B — > D\tv and B — > D\lv^ the former give rise to two 
independent universal functions, denoted by Ti(y) and r 2 (y) [Q. In the present work we 
shall focus on the first type of corrections and investigate these two form factors, Ti j2 (y), 
with QCD sum rules in HQET. 

The remainder of this paper is organized as follows. In Sec. |T| we review the formulae 
for the matrix elements of the weak currents including the structure of the Aqcd/^q 
corrections in the effective theory. The QCD sum rule analysis for the subleading Isgur- 
Wise functions related to the corrections from matching weak currents in QCD onto those 



in HQET currents is presented in Sec. |T|. Sec. [iy] is devoted to numerical results and our 
conclusions. 



II. THE HEAVY-QUARK EXPANSION AND THE SUBLEADING 
ISGUR-WISE FORM FACTORS 

The theoretical description of semileptonic decays involves the matrix elements of vector 
and axial vector currents (V 1 = c^b and A^ 1 = 07^756) between B mesons and excited 
D mesons. For the processes B — > D\iv and B — > D^lv, these matrix elements can be 
parameterized as 

(D 1 (v', e)| V» \B{v)) = f Vl + {f v ^ + f v ^) e*-v, (la) 
(D 1 (v',e)\A»\B(v)) = if A e^e* a v f3 v , 7 , (lb) 
(D* 2 (v', e)\A» \B(v)) = k Al e*^v a + {k A ^ + k M v'») e^v^ , (lc) 
(D* 2 (v',e)\V»\B(v)) = ik v e^6* aa v°Vf3v' 7 . (Id) 

The form factors fi and ki are functions of y — v ■ v', which can be expressed by a set of 
Isgur-Wise functions at each order in Aqcd / m c,b- I n the infinite mass limit a convenient way 
to evaluate hadronic matrix elements is by using the covariant trace formalism developed 



in Ref. [[03]. The heavy mesons can be described by spin wave-functions with well-defined 
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transformation properties under the Lorentz group and heavy-quark symmetry. The spin 
doublet parametrized by jf can be represented by a 4 x 4 Dirac-type matrix. For jf = l~/2 
and 3 + /2, the matrix representation are 



l + j> 



(2a) 
(2b) 



where P v , P*^ and P£, P*^ u are annihilation operators for members of the jf = 1 /2 and 
3 + /2 doublets with four- velocity v in HQET. The matrices H and F satisfy fH v = H v = 
-H v j>, j>FS = F£ = -FM, i^ 7 „ = 0, aadu M F^ = 0. 

Up to the order 0(l/m c &) the current cTb is represented in HQET by 



cTb = h$ ( r 



2m 



1 ^tT-i-i ^ 



2m& 



(3) 



where is the heavy quark field in the effective theory. Hadronic matrix elements of 
the leading order current between the states annihilated by the fields in H v and F", are 
written as 



h$Th® =T r Fr{v <T F^TH v }. 



(4) 



Here r is a universal Isgur-Wise function of y. 

At order Aqcd/^c,6 there are corrections originating from the matching of the b —>■ c 
flavor changing current onto those in the effective theory. For matrix elements between the 
states annihilated by the fields in F°, and H v , the new operators in Eq. (f|) at the order 
Aqcd / m c,b can be written as 



iD\ T hf = Tr {45 F V Z T H v } , 
r i~F)\ hf = Tr [S^l F°, Y H v ) . 



The most general decomposition of the form factor is 



S. 



(Q) 

aX 



' (Q) i (0) / i (Q) 



+ t 4 % ctA . 



(5) 



(6) 



The functions t, depend on y and have mass dimension one. Using the equation of motion 
for the heavy quark, iv ■ D = 0, and translation invariance, 
A'v' u ) htf T h$\ one can obtain the constraints H 



(7) 
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and relations between rj and r,- 

r^+rf^Ar, r< c > + r<» = -A' r , 

4 C) + = , r 4 (c) + ri b) = , (8) 

where mg + A and m<g + A' are masses of doublets jf = l _ /2 and 3 + /2 in the leading 
order. Furthermore, combining Eqs. (0) with ([8]) yields 

r 3 (c) =Fi (c) + r 2 (c) , 

ri c) = (y - 1) {t[ c) - r 2 (c) ) - (yA' - A) r. (9) 

These relations show that all corrections to the form factors coming from the matching of 
the weak currents in QCD onto those in the effective theory are expressible in terms of A r 
and A'r and two of the coefficient functions, say, and r 2 . 

Applying the above relations to the l/m c correction to the vector and axial- vector 
currents . and dropping the superscript c on r^, the corresponding matrix elements can 
be written as 

(D l (v', e)] ipr \B{v)) = i{ [ 4 (y + l)(A'y - A)r - 3 (y 2 - l)( Tl - r 2 ) ] e*" 

+{[ (4y - l)n + 5r 2 ] - [4(A'y - A)r + (y + 2)n 
+(3y + 2)T 2 ]^}e*- V ) , (10a) 



(Di(v', c)| ^ i^Tfe fc?> |S(«)> = ^[4(A'y - A)r - 3(y - l)(n - r 2 )] ^4V 7 , (10b) 

(£>*(?/, e) | i^ 7 " /#> = i ( n - r 2 ) £^ > ( 10c ) 

(£>*(?/, e)| ^ z^ 7 M 75 = (y - l)(n - r 2 ) e*^ v a + [ 2r x v» + (n + r 2 ) i/") ] 

x e^t;V. (lOd) 

The analogous formulae for the matrix elements of the 1 / my, correction to the currents can 
be written down in a similar way. 

There are also order AqcD/ m c,& corrections originating from terms in the HQET La- 
grangian of this order. They can be parametrized by additional functions of y [|]]. From 
now on we shall confine our attention to the determination of subleading Isgur-Wise func- 
tions, Ti(y) and Tz(y), associated with the matching of the vector and axial vector currents 
in full QCD onto those in the effective theory. 



III. SUM RULES FOR n AND t 2 

A basic element in the application of QCD sum rules to problems involving excited 
heavy mesons is to choose a set of appropriate interpolating currents in terms of quark 
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fields each of which creates (annihilate) an excited state of the heavy meson with definite 
quantum numbers j, P, je- The proper interpolating current J^pj" 3 for the state with 



arbitrary quantum number j, P, je in HQET was given in p0| . These currents have nice 



properties. They were proved to satisfy the following conditions 



\J&£ i (0)\f,P'Jt) = i fp^ PP ,5 jd > V a ^ , (11) 
f.f 

■J i 



J dt6(x -vt)U Pde (x) (12) 



in the itlq — > oo limit. Where r/ Ql '" Q j is the polarization tensor for the spin j state, v is the 
velocity of the heavy quark, gf = g al3 — v a v 13 is the transverse metric tensor, S denotes 
symmetrizing the indices and subtracting the trace terms separately in the sets (a% •••%■) 
and (Pi ■ ■ ■ Pj), f P j e and H P j t are a constant and a function of x respectively which depend 
only on P and j£. Because of equations (JTXj) and ([12]), the sum rules in HQET for decay 
amplitudes derived from a correlator containing such currents receive contribution only 
from one of the two states with the same spin-parity (j, P) but different je in the wlq — > oo. 
Starting from the calculations in the leading order, the decay amplitudes for finite rriq can 
be calculated unambiguously order by order in the I/tuq expansion in HQET. 



Following [^C§ the local interpolating current for creating pseudoscalar B meson is 
taken as 

J o,-,i/2 = \f\ hvlsq , (13) 

and the local interpolating currents for creating 1 + and 2 + (Di, D%) mesons in the doublet 
(Di, D%) are taken as 

-i) (v* - q , (14) 




Jt'.3/2 = \/i [iTW+lTW - J 9? ia2 ft) q , (15) 



'2, +l3 /2 - y \^ { ~y- {^ 1V ? 2 + ^v? 1 - 1 

where T> is the covariant derivative and 7^ = 7^ — ft) 1 * '. Note that, without the last term 
in the bracket in ( P^) the current would couple also to the 1 + state in the doublet (0 + , 1 + ) 
even in the limit of infinite tuq. 

The QCD sum rule analysis for the subleading form factors, Ti(y) and T2(y), proceeds 
along the same lines as that for the leading order Isgur-Wise function r(y). Here we shall 
briefly outline the the procedure for determining the Isgur-Wise function r(y) and refer for 
details to Ref. 0. The idea is to study the analytic properties of the three-point correlators 
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i 2 j Wz^) (0\T(r w/2 (x) J$ v '\0) 4 } _ >1/2 (z)) |0> = E{u,,u/,y) Ofa , (16a) 

z 2 / ^^e^-*-) (0|T(jg j3/2 (a;) J7^ y) (0) 4_, 1/2 (*)) |0> = S( W ,o/,v) ^ , (16b) 

where Jy- V ' v ^ = h(v')^ h(v) and J^ v,v "* = h(v') , y fl '-f 5 h(v) are leading order vector and 
axial vector currents, respectively. The variables k, k! denote residual "off-shell" momenta 
which are related to the momenta p of the heavy quark in the initial state and p' in the 
final state by k = p — mqv, k! = p' — rriQiv', respectively. For heavy quarks in bound 
states they are typically of order Aqqd and remain finite in the heavy quark limit. Cy,A 
are Lorentz structures associated with the vector and axial vector currents. 

The coefficient H(u;, to' , y) in (|TJj) is an analytic scalar function in the "off-shell energies" 
to — 2v ■ k and to' = 2v' ■ k' with discontinuities for positive values of these variables. It 
furthermore depends on the velocity transfer y = v ■ v which is fixed at its physical region 
for the process under consideration. By saturating the double dispersion integrals for 



the correlators in (16|) with physical intermediate states in HQET, one finds the hadronic 
representation of the correlator as following 

_ /_ if+,3/2 T (y) 

E hadro (uj, u', y) = -—= ' 2 . ' x : rv + higher resonances , (17) 

(2A_ i/a -u- ze)(2A +i3/2 - u' - xe) 

where fpj t are constants defined in flTTl), Ap,j t = mpj e — tuq. As the result of equation 
(|TT|), only one state with j p = 1 + or j p = 2 + contributes to (|T^) , the other resonance with 
the same quantum number j p and different ji does not contribute. 

Following the usual procedure of QCD sum rules and making double Borel transforma- 
tions in the variables u and a/, one obtains the sum rule for r as follows 

T(y) A 1/2 f+,3/2 e- 2 ( A -.V 2 m+A + ,3/ 2 /T 2 ) = J dudu 'pPert^ ^ y y-u / Tl -u> /T 2 

_<m_ a, f TiT 2 [(l + 2y)T 2 + 3T 2 + (Ay + 2)TiT 2 ] 
6 T 2 [ n '{ A%{T 2 + T 2 + 2yT 1 T 2 f 

(y 2 -l)(T 1 + T 2 )T?T 2 
12(T 2 + T 2 + 2yT 1 T 2 f 

where T\ and T 2 are Borel parameters and 

= ^{ 7 ^=j 3 ^- 1 l~^ 2 + (2y - 1)(-' 2 + 2-')] 

x 0(u) Q(u') &(2yvv> - v 2 - v' 2 ) . (19) 

The integration domain D is constrained by the conditions v < u c , v' < uo' c together with 
the functions in (|TD|) . 

Let us now turn to the derivation of the QCD sum rules for the subleading form factors, 
Ti(y) and r 2 (y), defined in (|6]). To this end, we consider the following three-point correlation 
functions 



| = K T (Ti, T 2 , u c , (J^ y) , (18) 
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* 2 j d*xd*ze< M ^ (0\T (jl + . 3/2 (x) h$0>fhP(O) Jj_, 1/a (z)) |0> = ^f/.y) , (20a) 

z 2 / rfVze^'- t2 » (0|T (Jr i+i3/2 (x) ^^^75^(0) j£- f i /a (*)) |0> = H^^f/.y) , (20b) 

* 2 / Wze^' (0|T(j 2 ^ 3/2 (x) ^z^/i^O) j£_, 1/a (*)) |0> = E%f(u>,u/, V ) , (20c) 

z 2 / Me*'' 1 ^ (0|T(jg )3/2 (x) ^^^75^(0) 4,-,i/ 2 (*)) \0) = E&(u>,u/,v) ■ (20d) 



By saturating the double dispersion integral for the three-point functions in Q2"0[ ) with 
hadron states and using ([II]) and ([D]), one can isolate the contributions from the double 
pole at uj — 2A_ 1 / 2 , uj' = 2A +i3 / 2 : 



/_ i/ + ,3/2[4(yA + , 3 /2 - A_, 1/8 )r(y)£& T + r^y)^ + r 2 (y)C 



[IV 

1Vt 2 1 



~iv(w,o;,y) = ? j-= — -= + •••, 21a 

(2A_ ) i/ 2 -uj - ze)(2A +i3/2 - uj' - it) 

„, VI , ^ /- 1 i/+,8/a[4(yA + ^ /a -A_, 1/a )r(y)-3(y-l)(r 1 (y)-7i,(y))]/:^ 

a$ A {w,w,y) = 2 - —= — tt ; ; 7 + •■•, 21b 

(2A_,i/ 2 -uj- ze)(2A +i3/2 - uj' - le) 



wp, , v /-,§/+,3/2[ri(y)-r 2 (y)]£^ 

i$r(u,w ,y) = j-j ? — 77TT : -- + ■■■ , (21c) 



,3/2 [ 

(2A_,i /a - w - ie)(2A +j3/2 - u/ - ie) 
/-,^/ + ,3/2[ri(y)£^+r 2 (y)£ffl 
(2A_ il/2 -cu-ie)(2A +i3/2 



^,u',y)= 2 5 "-'-;; w ;^;' + . (2 i d ) 



where the Lorentz structure £ M ^'s and £ MQ/3, s are collected in the Appendix. 

On the other hand, it turns out that the calculation of these correlators with the operator 
product expansion (OPE) in HQET yields the following general structure 

Ely(uj, uj', y) = E f Clvr + S n^iy T1 + z^Cwn > ( 2 2a) 

E^(uj,uj',y) = [E, - 3(y - l)(H n - H T2 )]^ , (22b) 

E&{u,c/ t y) = (R n -E n )/%* , (22c) 

SgV a/, y) = S T1 £^ + ^ . (22d) 

where the coefficient functions 5 T1 , S r2 and Hf are scalar analytic functions in the off-shell 
energies uj and uj'. 

Comparing (21) with (22) one can see that they are compatible. Therefore, we can 
calculate the scalar functions Sj(u;,u/, y) with the QCD sum rules in HQET. From (|22"D 
and (PH) one can see that the sum rules for H Tl , S T2 and yield sum rules for Ti(y), r 2 (y) 
and (yA +)3 / 2 — A_ 1 / 2 )r(?/), respectively. In the theoretical calculation, for simplicity, the 
residual momentum k is chosen to be parallel to v such that = (k ■ v)v^ (and similar 
for k'). The theoretical expression for the correlator in HQET consists of a perturbative 
part and contributions from vacuum condensations. Confining us to the leading order of 
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perturbation and the operators with dimension D < 5 in OPE, the relevant Feynman dia- 
grams are shown in Fig 1. We shall focus, at first, on the coefficient functions H n (u;, a/, y) 
and S T2 (c<j, lo\ y) to construct the sum rules for the subleading form factors Ti(y) and r 2 (y). 

The spectral densities in the double dispersion integral for the perturbative diagram 
depicted in Fig. 1(a) turn out to be 



2 8 vr 2 \Vy 2 - V y + i 

-(2y 2 -2y + l)v' 3 



2% 2 VVy 



y+l 



5z/ 3 + (12y - 3)z/V - (6y 2 - 6y + 3)zV 2 
0(v) 0(i/) Q(2yuu' -v 2 - J 2 ) , 

(Ay - l)z/ 3 - (9y 2 - Qy)v 2 v' + (Ay 3 - 8y 2 



+2y-l)vv' 2 + (2y 3 -y 2 + 2y)v 



'3 



Q(v)Q(v')Q(2yvv' -v 2 -v 



'2\ 



(23) 



(24) 



The non-perturbative power corrections to the correlators are computed from the diagrams 
involving the quark and gluon condensates in Fig. l(b)-(c) in the Fock-Schwinger gauge 
x fl A tl (x) = 0. We find that the only non- vanishing contribution is the gluon condensate. 
After adding the non-perturbative part and making the double Borel transformations one 
obtains the sum rules for T\(y) and r 2 (y) as follows 



Ti(y)/- l i/ 2 / + ,s/ 2 e- 2 ( A -.v 2 m+A + , 3/2 /T 2) 



D 



dvdvY T ert (v, v\ y)e- u/Tl - u ' /T2 



(—GG)H T1 (T X ,7 2 ) 



71 



r a (y) f-,1/2 / + ,3/ 2 e- 2 ( A -.v 2 /^+ A + ,3/ 2 /^) = f dvdv'f? T f(v, p>, y)e^/ T ^'^ 



D 



■(—GG)H n (T h T 2 ) 

7T 



(25) 



(26) 



where 



1 



H T1 (T U T 2 ) = -T?T 2 



Tl + AyT^T 2 + (Sy + U)T?T 2 + (Ay 2 + 12y + 1Q)T 2 T? 



+(16y + 3)T!T 2 4 + ATI 



/(T 2 + 2yT 1 T 2 + T 2 ) 4 



(27) 



H T2 (T U T 2 ) = -T 2 T 2 



(2y + 1)7* - (8y 2 - 8y)T?T 2 - (28y - 9)7?7 2 + -2Ay 2 7 3 7? 



-(8y 2 + 6y- 11)T 2 T 2 4 + 8y7{7^ + 3T, 



j(7l + 2yT x T 2 + T| 



(28) 



The integration domain D is restricted to the area in v < uj c , v' < oj' c . 

We have checked that the sum rules for T\(y) and T 2 (y) derived from Eqs. fl22aj) and 
( |22d|) are the same and they are also consistant with the sum rules derived for t\ — r 2 from 
Eqs. O and ( ggg) . 

Furthermore, from the coefficient function Hf in ( |22a| ) and (|22b| ) one finds the same 



sum rule for the combination (yA_ 



,3/2 



A_ ; i/2)t(2/) as follows 
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(2/A + ,3/ 2 - A_ |1/2 )r(y) /_, 1/2 / +>3/2 e- 2 ( A -^ 2 /^+A + , 3/2 /T 2 ) = 

~2 (^dT^ 1- Q^-[) K r{Ti,T2,uJ c ,u' c ;y) , (29) 



where X T in this equation is identical to the function K T in (|18D . Therefore, ( P9|) is 



consistant with the sum rule (18) in the leading order, 

The above consistancy checks confirm that the correlators have the forms in (22) and 
that our method is consistant with the general analysis of Ref. Q described in Sec. |I[ 

In obtaining the sum rules fl2"5|), (p6|) and (^) the quark-hadron duality has been 
assumed. In doing this, the contribution from higher hadronic states is simulated by the 
perturbative part above some threshold energy. In the QCD sum rule analysis for Isgur- 
Wise functions for B semileptonic decays into ground state D mesons, it is argued by 
the authors of |2||||7j that the perturbative and the hadronic spectral densities can not be 
locally dual to each other, the necessary way to restore duality is to integrate the spectral 
densities over the "off-diagonal" variable = [y — v')/2, keeping the "diagonal" variable 
v + = (z/+z/)/2 fixed. It is in z/ + that the quark-hadron duality is assumed for the integrated 
spectral densities. In doing this, for simplicity, the two Borel parameters are taken to be 
equal: T\ = T 2 = 2T. We shall use the same prescription here. 

The functions in (^) and (|24|) imply that in terms of v + and z/_ the double dis- 
continuities of the corrrelator are confined to the region — y/y 2 — 1/(1 + y) v + < v_ < 
\/y' 2 — 1/(1 + y) u + and v + > 0. According to our prescription an isosceles triangle with 
the base v + = v c is retained in the integation domain of the perturbative term in the sum 
rule. 

In view of the asymmetry of the problem at hand with respect to the initial and final 
states one may attempt to use an asymmetric triangle in the perturbative integral. How- 
ever, in that case the factor (y 2 — 1) 3//2 in the denominator of (|2~3"D and (^4j) is not canceled 
after the integration so that the Isgur-Wise function or it's derivative will be divergent 
at y — 1. Similar situation occurs for the sum rule of the Isgur-Wise functions for the 
tansition between ground states if a different domain is taken in the perturbative integal 

0- 

Putting everything together one obtains the final expressions for the QCD sum rules 



1 ^ GG ) y + 9 T , (30) 



3x2 5N tt (y+l) 3 

f\t f ~( A i+ A +3)/r 1 5y-l 4 /T 

2 K yjj-,-j+,- 16tt 2 V + l 4 Jo 



1 -(— GG)-^ — |- T . (31) 



3 x 2 5 v 7T ' (y 
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We end this section by noting that the QCD 0(a s ) corrections have not been included in 
the sum rule calculations. However, the Isgur-Wise function obtained from the QCD sum 
rule actually is the ratio of the three-point correlator to the two-point correlator. While 
both of these correlators subject to large perturbative QCD corrections, it is expected that 
their ratio is not much affected by these corrections because of cancelation. This has been 
proved to be true in the analysis for B semileptonic decay to ground state heavy mesons 
0- 



IV. NUMERICAL ANALYSIS AND CONCLUSION 



In order to obtain the values for Ti(y) and T 2 (y) from Eqs. (|30|) and (|31|) in the numer- 
ical evaluation we need to use the hadronic parameters A's and /'s of the corresponding 
interpolating currents as input. The QCD sum rules derived from the two-point correla- 
tor has been applied to determine A and /. A—,1/2 and f-,1/2 can be obtained from the 
results in [0] as A_ 1/2 = 0.5 GeV and /_ 1/2 ~ 0.24 GeV 3 / 2 at the order a s = 0. No- 
tice that the coupling constant f-,1/2 defined in the present work is a factor l/y/2 smaller 
than that defined in 0. Determination of A +j3 / 2 and /+ j3 / 2 by QCD sum rules gave the 
results: A +j3/2 = 0.95 GeV and / +)3/2 = 0.19 GeV 5 / 2 at the order a s = 0,0]. For 
the QCD parameters entering the theoretical expressions, we take the standard values: 
(a s GG) = (0.04) GeV 4 . 

Imposing usual criterium that both higher-order power corrections and the contribution 
of the continuum should not be very large, we find an acceptable stability window for 
the threshold parameter in the range u c = 2.0 — 2.6 GeV, in which the results do not 
appreciably depend on the Borel parameter in the range T = 0.7 — 1.1 GeV. The range of 
Borel parameter here overlaps with that of the corresponding two-point sum rules |7|j20| . 



The values of the form factors T\(y) and r 2 {y) at zeor recoil as functions of the Borel 
parameter are shown in Fig. 2(a) and 3(a), for three different values of the continuum 
threshold uj c . The numerical results for T\{y) and T2(y) are shown in Fig. 2(b) and 3(b), 
where the curves refer to three different values of uj c and T is fixed at T = 0.9 GeV. 

The numerical analysis shows that T\{y) and T 2 (y) are slowly varying functions in the 
allowed kinematic range for B — ► D\lv and B — * D\iv decays. The resulting curves for 
T\{y) and t 2 (?/) may be well parameterized by the linear approximations 



n(y) = n(l) (1 - p 2 Tl (y - 1)) , n(l) = -0.4 ± 0.1 , p 2 i = 1.4±0.2, (32) 
r 2 (y) = r 2 (l) (1 - p 2 T2 (y - 1)) , r 2 (l) = 0.28 ± 0.05 , p 2 T = 0.5 ± 0.1 . (33) 



The errors here reflect only the uncertainty due to u c and T. They do not contain other 
errors in the QCD sum rule approach. 
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In conclusion, we have presented a QCD sum rule analysis of the subleading Isgur-Wise 
functions Ti(y) and T2(y), appearing in the heavy quark expansion of the transition matrix 
elements between heavy mesons due to matching of the weak currents in QCD onto those 
in the effective theory at the order 1/mq. Our approach is in accordance with the general 
relations obtained from analysis based on HQET in |4j]. 
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APPENDIX: 



We list here the lorentz structures used in the paper. 

C^ VT = ^=[-(y + l)gr + v'^}, (Al) 

Wv n = ^{s(y 2 - i)gT - [(% - iK - (y + 2)v»w} , (A2) 

C& n = ^{- 3(y 2 - i)gT - [W - (3y + 2)v'w) , (A3) 



CZ = ^ af3 v a v'p , (A4) 

Q$ = + ^ P ^>X , (A5) 

3S£ = (y - l){i(^f + flfV) " + (2^ + v'n{v?v? - i(l - y 2 )^} , (A6) 

= -(i/ - i){^(sT«f + sfV) - \af<} + - 5(1 - y 2 )^} , (A7) 

where g? = g alS — v' a v'P and t>f = v a — yv' a . It is easy to see that these Lorentz structures 
satisfy 



°u L '1Vt — V v L -\Vt x — U v L -\Vt 2 


— V u L 'lA 


= 


(A8) 


ga/3J~-2V — QaP^A-n ~ 


9a/3^2Ar 2 


= 


(A9) 


V a L '2V — V a L -2Ar 1 ~ 


„' r^P . 

u u'~'2At2 


= 0, 


(A10) 


1 r iia(3 _ 1 rl iaP _ 
U /3 J ~'2V — u P l '2At 1 ~ 


- V P L '2At2 


= . 


(AH) 



The appearence of the Lorentz structures satisfying these relations is the result of the 
following equations 

< Jt+,3/2 = < Jf +t3f2 = , (A12) 

tq/3 n ja(3 jfla 1 a -I o\ 

9a/3 J 2,+,3/2 ~ U ' J 2,+,3/2 ~~ J 2,+,3/2 ' \ Al6 J 

satisfied by the interpolating currents in the correlators. 
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Figure Captions 



Fig. 1. Feynman diagrams contributing to the sum rules for the Isgur- 
Wise form factors in the coordinate gauge. The gray square corresponds 
to the insertion of the l/m c current. 

Fig. 2. Numerical evaluation for the sum rule (]30|): (a) dependence 
of Ti(l) on the Borel parameter T for defferent values of the continuum 
threshold uj c ; (b) Results for the Isgur-Wise form factor T\{y) with T = 
0.9 GeV. 

Fig. 3. Numerical evaluation for the sum rule (pT|) : (a) dependence of 
T2(l) on the Borel parameter T for defferent values of the continuum 
threshold u> c ; (b) Results for the Isgur-Wise form factor Tziy) with T = 
0.9 GeV. 
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